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ABSTRACT 


A boundary value problem of two elastic bodies in 
contact is considered. The bodies are finite circular cyl- 
inders of different dimensions and material constants and 
are isotropic and homogeneous. They are forced into contact 
across the plane faces such that the resulting stresses and 
displacements are axisymmetric . The solution utilizes 
Love's stress function to generate a family of biorthogonal 
eigenfunctions for each cylinder. The interrelated Fourier 
coefficients are expressed implicitly by an infinite system 
of linear algebraic equations. By truncation, an explicit 
solution of the Fourier coefficients is obtained. 

Two example problems are solved* first, the two 
cylinders are placed in frictionless contact; second, the 
two cylinders are placed in bonded contact. In each problem 
the other face of each cylinder undergoes a constant displac 
ment with zero shear tractions. Selected numerical results 
are presented in graphical form. 
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CHAPTER I 


THE PROBLEM 

In this thesis, a method for solving the contact 
problem of two finite circular cylinders is presented. 

The cylinders are forced into contact across the plane 
faces such that the resulting stresses and displacements 
are axisymmetric. They are compressed between two rigid, 
frictionless faces as shown in figure 1. 

Both cylinders are assumed to be homogeneous, 
isotropic, and elastic. Each cylinder may have its own 
radius, length. Young's modulus, and Poisson's ratio. 

The curved surfaces will be assumed to be free of trac- 
tions. 

Two example problems will be solved: first, the 

two cylinders are placed in frictionless contact with each 
other; second, the two cylinders are placed in bonded 
contact. 

The solution will be based on Power [1] and Power 
and Childs [2]. Other solutions for the finite cylinder 
do exist but are not satisfactory because of numerical 
properties or approximate schemes for meeting the bound- 
ary conditions. 




CHAPTER II 


THE FINITE CYLINDER 

The solution for a finite cylinder problem* as 
formulated by Power M. is based on a Love function [3] 
in the form of a series of Bessel functions times expo- 
nential functions plus an improper Fourier integral. 

The Love function was taken to be the sum of three solu- 
tions 


X = X° + X'+ X" (2-1) 

where X° is the zero eigenvalue solution? X ’ is the self- 
equilibrated solution associated with the complex eigen- 

p 

values; and X is the particular solution associated with 
body forces, which will be neglected. 

The governing equation for the homogeneous solutions 
X ° and X 1 is 


V 4 X = 0 (2.2) 

Each of the separate solutions satisfies the boundary 
conditions on the curved surface. 

In terms of the Love function, X, the stresses and 
displacements that are not zero are given by [3] s 
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— 

[(1-WV ! X- X, ZI ] 

s- j t r 

(2.3d) 

E u 


-*. rI 

( 2 . 3e ) 

1 + V 
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(2.3f) 
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The 

zero 

eigenvalue solution of the 

biharmonic 


equation V 4 X =0 is found by separation of variables. 
A sufficient portion is 


X° = Az 3 + Bz 2 + Ar 2 z (2.4) 

1-2 V 


For a cylinder of radius R, the stresses and displacements 


v 0 

obtained from A are 


cr° = c r;= r°= 0 

(2.5a) 

^0 _ 6(1 + V) A 

^ " 1 -2V A 

(2.5b) 

E u° 6 V r A 

(2.5c) 

(1+V)R 1-2 V R 

E w° _ 6 z A , „ B 

(1-f l/)R ~ 1-21/ R A + R 

(2.5d) 
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It is more appropriate to call this solution the 
uniform axial stress solution, as can be seen by equation 
(2.5b). Equation (2.5d) shows that this solution produces 
a net axial displacement consisting of two terms. The first 
term is proportional to z. The second term is the net axial 
displacement at z = 0 and is considered as a rigid body 
displacement . 

The self-equilibrated solution, X 1 , for the finite 
cylinder is based on Little and Childs solution of the 
semi-infinite cylinder. This solution, containing only the 
eigenvalues in the right half of the complex plane, was 
extended by Power [l] to solve the finite cylinder by using 
the eigenvalues in the left half of the complex plane and 
the axial coordinate z- L . This solution for a finite cyl- 
inder of radius R, length L, and with self-equilibrated 
tractions on both ends is 
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The dimensionless eigenvalues, < Y ] , are defined as 

^ d A 


(2.7) 


where the £jR are the roots of the characteristic equation 

UR) 2 J 0 2 UR) + [(£ R ) 2 - 2( 1-V)] J 2 (*R) = 0 (2.8) 

The eigenfunctions, <p (Yj , r/R), are also dimensionless and 
are given in the Appendix. 

The specified boundary values of X^are defined as 


(r) = cTj( r, o) = 


cr b - cr ( 


x = o 


(2.9a) 


h, ( r) = f'( r, o) = T 


2—0 


(2.9b) 


q. ( r) = — m* (r, o) - — l"u b - u° 
j 1 i + v J i + v L 


2 = 0 


(2.9c) 


K, ( r) 


i ^ > i/ E» 


I + v 


w '( r, o) ~ j-pg- w b - v \/ 0 


[« 


2=0 


( 2 . 9d ) 


fUr) = <r 1 ( r, l) 


2 = L 


( 2 . 9e ) 


h Jr) = r 1 ( r j l ) = r b 


Z. - L 


( 2 . 9f ) 
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g 2 (r) 

L 

E. 

11 

V? 

V* 

3 

£ 
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i 4* U 

1 + V 

2 -u 

K z (r) 
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( r, l) = 

E 

w b - w ° 


(2.9*0 


l +• V 

I + V 

2 = L 


The Fourier coefficients are therefore obtained from 
the integrals 

ft 

a ij = n crj ) f [Wi(tji^)£Cr;+W l (i5 J ^)h l (r) + 

o 

+ W,(r4)K,Cr)]lf 

P 

Q ij = ^)J E W ' (t* W ^'S) ' 

0 

+v^(r i ^)K i (r>] £ £ 

where 

N(7j ) = 2(l-V)[z(l-V)J„(rj)J,(7j) - Z^JoCTj) 

Explicit forms of the W functions are given in the 
Appendix. 

Only four of the functions (2.9) may be arbitrarily 
specified as self-equilibrated boundary conditions. The 
other four are obtained from equations (2.6) evaluated at 
the appropriate boundary. Equations (2.10) may then be 
written in the form 


f,C r) 


(2,10a) 


feCD 


(2,10b) 


Z ^Z )] 

( 2 . 11 ) 
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K = 1 
GO 


K- I 
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^ aj “ ^ Z j 2D ^ZjK^2.K^"2D ^2.j K ^ I ! 
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(2.12a) 


(2.12b) 


where Gj_ j and G 2 j represent that part of the integral con- 
taining the specified boundary conditions on X\ 

It should be noted that, in order to avoid a null 
self-equilibrated solution, the specified boundary con- 
ditions must contain at least one of the functions (2.9) 
where the uniform axial stress solution does differ from 
the boundary condition. 



CHAPTER III 


TWO CYLINDERS IN CONTACT 

The solution for the two cylinder contact problem 
will be found by using the finite cylinder solution appli- 
cable to each of the two cylinders. The equations will then 
be solved simultaneously to yield the Fourier coefficients 
in accordance with the boundary and continuity conditions. 
The geometric dimensions of the two cylinders are indicated 
in figure 1. For the lower cylinder, its eigenvalues and 
Fourier coefficients will be denoted by y. and a n and a 2j , 
respectively. Similarly, X, and a 3j and a 4j will be used 
for the upper cylinder. Subscripts ( ) L and ( ) y will 
denote other values for the lower and upper cylinders, 
respectively. 

The Uniform Axial Stress Solution . 

Equations (2.5) from the previous chapter will be 
used. It is convenient to write them in this forms 


^zt- = „ j 0~zu - „ A u (3»la) 

Al $ Al A l 


A 





(3»lb) 


Jk^k = - z + = z — +- — 

A l A u j A l Al Al 


(3-lc) 
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The Total Solution * 

The total solution will be the sum of the self- 
equilibrated solution and the uniform axial stress solu 
tion. For the lower cylinder, it is 
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oo 
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oo 

E 
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oo 
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Similar equations can be written for the upper cylinder, 
An example is 


Al 


j*i 


Au 

A l 


The Fourier Coefficients . 

The following functions will be defined for speci- 
fying the boundary values of X' for the lower cylinders 


P.(r) 


0 * ' ( r, o) 


q,(r) 


E l 


0 + VJ 


?Ar) 


(T^Cr.i) 


Sz< r ) 


E-l 


C I + vJ 


(3.2a) 
(3.2b) 
(3-2c) 
(3- 2d) 
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Ut( r, i) 
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h i(r) = ri < r 3 o) 5 h (r) = V^C^i) ( 3 . 4 c) 

K,CD = w; C r, o) , K z (r) = ^ ^y^rJ) (3-4d) 

Similarly, these will be defined for the upper cylinder 

f 3 oo = <Cc 00 ; , f/rj = O^cr, h) (3-5a) 

9 3 (rj = <rf%U>-,o) , g„(r) = <^u;,<r,h; (3-5*» 

h 5 (r) = ru(r, O) , h,(rj = rj (r, h) (3-5=) 

K 3 ( r) = ( T R 1) '*t (r ' 0) > l< v (r) = of% w u'( r - h) (3.5<J) 
The <f> and W functions exhibit an inner bi orthogonal! ty 
in the form 

i 

nT*T /[ w n(t, r)0 n (r K ,r) + W nt2 (Tj,r) $ n+ p K ,r)} r dr 

J O 

= y 3jk * n = 2. (3.6a) 

wfe 1 W»- 4 ) 1 *~M + 4 ) «wM)] 4 f 

= x 4 .« ■ 


n - 1 , 2 . 


(3.6b) 
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Using this property, the Fourier coefficients are obtained 
from the integrals : 

1 

a 'j = myr/ [w,(Yj,r) fi(r) -t- W 2 (rj,r) h,(r) 

' J ' 0 

+W 3 ( 7 j,r) ^.(r) + W 4 (%-, r) K ,(r)] rdr (3. 7 a) 

1 

a 2J = N(^/[~ w i(Yj,*')f l (i r )+ W z (-r/.n 

-W 3 (Yj,r) g 2 (r) + W lf (-Y|,r) K z (r)] rdr (3.7b) 

V 

a 3j = i?fe/[w.K£)6(r) + W 8 (^)h 9 (rJ 

J O 

+W 3 (>j,^) 9 3 (r) K s (r)] rdr <3.70 

Q - , tj = n ( Aj) / F^ r ) + W 2 h»(rl 

-W 3 (^4) j^rJ+W^^r)] rdr <3.7d) 

It can be shown that when the specified boundary con- 
ditions pertaining to a two cylinder contact problem and the 
unspecified conditions are substituted into equations (3.7) 
and then integrated, the Fourier coefficients can be expressed 
in the form of the following matrix equations: 
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(3-8a) 


(3.8b) 


(3.8c) 
( 3« 8d) 


where G m j , the forcing function vectors, are that part of 
the integral (3.7) which contains the specified boundary 
conditions on X’. The matrices and Cnjif for n = 1 

through 4 represent the part of the integral ( 3 . 7 ) contain- 
ing the unspecified conditions on X The unspecified conti- 
nuity conditions are represented by the matrices £5^, 

c 5 jk> C6jk‘ Unspecified conditions are based on the 
solution for the stresses and displacements evaluated at 
the boundary. Generally, the infinite systems of equations 
(3.8) can be satisfactorily approximated by truncation. 
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As stated by Power [ 1 ] , this would be true if the matrices 
are strongly diagonal, which would result in the lower 
order coefficients being only slightly influenced by the 
higher order ones. The matrices B n j ^ and C n jj c for n=l 
through 4 does exhibit this characteristic. However, the 
matrices B 6 jk, and G 5 contain eigenfunctions 

from one cylinder and biorthogonal functions from the other 
cylinder. It is believed that they are somewhat strongly 
diagonal depending on the configuration of the cylinders 
and material constants- A positive proof is beyond the 
scope of this thesis but it will be assumed that all ma- 
trices can be truncated to yield results that are suffi- 
ciently accurate. 



CHAPTER IV 


EXAMPLE PROBLEMS 

Two particular problems will be considered! 

Problem 1. The cylinders are in frictionless contact 
with each other. This implies the continuity of normal 
stresses and normal displacements. In addition, the radial 
shear stress will be zero at the interface boundaries. 

Problem 2. The cylinders are bonded together. This 
implies that the normal displacements, normal stresses, 
radial displacements, and radial shearing stresses are con- 
tinuous through the contact interface. 

As a prelude, a short discussion of the indentation 
of an elastic half space by a rigid, flat-ended cylindrical 
punch will be made. The solution is concise and well known. 
It has many similarities to the two cylinder contact pro- 
blem. For a rigid cylinder of radius 7 ) indenting an elastic 
half space, the contact pressure is given by Sneddon [6] as: 


?i_ m - V 

p 2(^ 2 - r 2 F 


(4.1) 


where p is the average pressure based on the total load, F 



P 


(4.2) 
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Equation (4.1) shows that the contact pressure increases 
from §p at the center of the punch to infinity at the edge 
of the punch. 

It is expected that a stress singularity of this 
nature will exist for the two cylinder contact problem. It 
would occur at the interface for the value r = v • 

Boundary Conditions at the Rigid Faces . 

It is convenient to first consider the boundary con- 
ditions which are applicable to both problems. It will be 
recalled that the rigid compression faces are assumed to 
contact the cylinder ends in a frictionless manner. More- 
over, the bottom face remains stationary and the top face 
is displaced downward an amount equal to 5. Four boundary 
conditions for the two cylinders ares 

t L (r , o) = w L (r> o; = r u (r, h; = o ( 4 . 3 a) 

W g (r, h) - - S (4.3b) 

The related boundary conditions on X ' are then; 

h, (r ) = h 4 (r) = k,(r) = k 4 (r) = 0 (4.4) 

Substitution of these specified conditions and the appro- 
priate unspecified conditions into equations (3* 7a) and 
(3*7d) results in 
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11 

G,^- = O 

(4,5a) 


&*ytc “ ^ ^y'K 

(4,5b) 

^ l|K — 

- / 2 S y , e‘ rK ' e 

(4.5c) 


- ‘/ 2 s,-K 

(4. 5 d) 

These results are 

substituted into equations 

( 3 . 8 a) and 

( 3 . 8 d) to yield 



IK. 

= — e "”*** ^ q. 

(4.6a) 

O'** 

= - e~** h/ V a 3K 

(4.6b) 


It is apparent that these last two equations can be substi- 
tuted directly into equations (3* 8b) and (3* 8c). The 
results ares 



G 2 / + 

^K^3K 

(4.?a) 

3^‘K ~ 

s + 

^ 2 * 

(4.7b) 


where 
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- 

+ 

^ -Tv JL 

C <o K 

(4.8a) 

3 jk, = 

S/k. " 

B 3,K + 

C 3/ k 

(4.8b) 

«/■* = 

e>s,-K 

_ ^ 5 ft* 

g -^Vth/'Tj 

(4.8c) 

67 k. 




(4. 8 d ) 


Substitution of equation (4.7b) into (4.7a) results in this 
explicit matrix solution for a 3k i 

a 2 = [d, - D 3 D- 3 'D 6 ] [g 2 + D s D*3 G 3 ] (4.9) 

The other sets of Fourier coefficients can now be readily 
found by using equations (4.7b) and (4.6). 

The Constants for the Uniform Axial Stress Solution . 

The total solution for the axial displacements at 
the boundaries will first be given. Equations (4.6) and 
the appropriate value for z will be substituted into equa- 
tions ( 3 * 2 d) and ( 3 * 13 ) to give 


E. t w L (r ? o) 
A t (i+V u ) 


A u l+V L 


(4.10a) 


£-i Wj, ( r, £) __ r q . 1 (p (' r ■ f) - — ! — j. — L... (4.10b) 

A l ( I + V L ) “ 4^ [ u, j u *i J J \+v L Au i+K 

j = i 
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EijW u (r ,o) _ f>L 
AJI+V U )^““ ^L a »j 

E u w a (r,h) 

A l (i + Ui)*J A l 


\ h > 

1+ V u ^ A u l + V u 


+ 


I 


'n 


A t 


_J i_ 

H-Vu *| 


(4.10c) 

(4.10d) 


The boundary condition that w L (r, 0) = 0 is substituted into 
equation (4.10a) to yields 


B l 


= 0 


(4.11) 


The force equilibrium condition requires that the following 
be true 


/ tf ZL r dr = / CT zu r dr 
o o 


(4.12) 


Substitution of equations (3«2a) and (3*3) into (4.12) 
leads to 



(4.13) 


The continuity of displacements at the interface can be 
ensured by 

w L (r» l) - w u (r, 0) (4.14) 

This equation is valid for all values of r<r>. In applications, 
a finite number of the eigenfunctions will be used and it is 
then desirable to use the integral of this equation 
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p P 

Jw t (r, i) r dr = f w g (r, 0) r dr (4.15) 

o o 

where 0 <p<r). Substitution of equations (4.10b) and (4.10c) 
into (4.15) leads to 


where 



(4.16) 


C [ Q|J e Tj/+ - X C-rjPj-f^.CTjP) 

j ~> ' 

+■ P x, O'} p) - o + v„) ^ j t, ( »J ) |t-[a v + a %j e'* j h/r *] 


rz vuP 
Laj 





(4.17) 


The constant, A L , can now be determined. The bound- 
ary condition that w y (r, h) = - 5 and equations (4.13) and 
(4.16) will be substituted into equation (4.10d) to yield 


a l = 


E l 6 


h E l 
t? 2 Eu 


+ £ 


(4.18) 


The zero-eigenvalue solution is therefore found to be 


<r ° 

£*=- = - > 

** u t 

CThu _ _ _L_ 
A t " ^ 

(4.19a) 

= r V u 

A u * 

EuuC __ r v u 
Au ' v z 

(4.19b) 

A u > 

E-uWu _ Z EjjL.4-* 

A L " ->7* + E.^ •*) 

(4.19c) 
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Problem 1. Frictionless Contact # 

Consider a cylinder of radius 1 as being placed into 
frictionless contact with an upper cylinder whose radius, 7? * 
is equal to or less than 1. For the lower cylinder, the 
following boundary conditions will be used: 


cr, L (r,i) = cr au (r, 0 ), 

0^ 

(4.20a) 

3-0 ~ O 9 

7 ^ r ^ l 

(4.20b) 

r L (r,^; = o 


(4.20c) 


For the upper cylinder, these boundary conditions will be 
used: 


W„(r, o) = W L (r,j) 


(4.21a) 


T u (r, o) = o 


(4.21b) 


The self-equilibrated boundary conditions on X ' are thens 

fa = t K - e ~*' h/1 ] <R {h> 4) - - ] , 

o ^ r £ (4.22a) 


f*(r) = i 


f 


*7 ^ r ^ 1 


(4.22b) 
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h 2 Cr) 

= h 3 (r) = O 




(4.22c) 

K,fr) 

_ 1 E. 

“ T E- 

J,K eV+a! 

J <a^ r ) 

(4, 22 d ) 

The unspecified conditions are 
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(4.23b) 

f 3 (r) 

00 
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[a* - a 4 -e^ h/ \ 
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*) 

(4.23c) 

9.CO 

QO 

=21 

j=i 

[a« + 


^JJ 

3 

(4.23d) 


Substitution of these conditions into the Fourier integrals 
( 3 . 7 b) and ( 3 - 7 c) result in the following forcing functions 
and matrices: 

rjrdr - 

" J 9 o o 

" Q 


J W, r) r drl (4.24a) 
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N CX ) 


1 


(4.24b) 
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(4.24e) 


(4.24f) 


(4. 24g) 


(4.24h) 


(4.241) 


(4.24j) 
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Problem 2, Bonded Contact . 


Consider two cylinders as being placed into contact 
as in problem 1 except in this case they are bonded together. 
For the lower cylinder, these boundary conditions will be 


used: CT IL (r,A) = (Tju (r, 0) } 0 « r s 

(4.25a) 

t L (r , A) = t u (r,o) , 

(4.25b) 

CT il (r, i) = o , «) ^ r ^ i 

(4.25c) 

r L cr, = o , ss r s i 

(4.25d) 

For the upper cylinder, these boundary conditions 
used s 

will be 

Wu (r, o) = w c (r, X) 

(4.26a) 

u u cr, o) =■ u L ( r, x; 

(4.26b) 

The self-equilibrated boundary conditions on X ' in 
case are: 

this 

_ s e " :vh/ ’ ! ] (\4) ~ B 

i~' 

3- 1 

T j > 

o s r « ^ 

(4.27a) 

f g (r) = I ; ^ r ^ i 

(4.27b) 
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=2[ Q 3V + 


j o « r « 
r * 

(4.27c) 

h 2 (r) = o } * r i 


(4.2?d) 

9* (r) - e_ l i +v u 2 C a, / ( 

r l 

e 

<p 3 o;-, o 

<r* v v L r rvu e-l > 

r _ i”l J^iiL 

(i*.2?e) 

C + Lv l E-u^ 

f U E-u 

«. (r > = i l° L .'tv;2 [<v 

e + a 2 - 




(4. 2?f ) 

The unspecified conditions are 

} 


9i (r) =2[ a 'i e_v - 
r 1 

( v 4 , r ) 

(4.28a) 

°o _ 

Kjf) + a 2£ 

r< 

<P»OJ, r ) 

(4.28b) 

f, = 2[ a *- i V*'"*' 

i-< 

«(»4) 

(4,28c) 

h 3 (r) = 2 [ a *< + S e " A ' h/ ' n 
i~' 


(4.28d) 


Substitution of these conditions into the Fourier integrals 
as before result in the following forcing functions and 
matrices i 
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/rj I 

S' = Sik) &/ w - (r /» r)rdr -fW'(y,r)rdr] 

* O O 



- Vl 

E-u fVu i 

Lk ! 

_ ,1 r!./ /-> ^ r* dr 

q h 

NCA^O + VJ 

E-L Lv u E 

L u ^ 

J J vv 3 p’»-ny r 

o 1 

I 

n-k 

= N <r,)/[ W » 

* O 

(r j t r) <p 3 (x, r) +■ r ) <pj,.(X, r) 

rdr 


e“ Y,<£ rr 

=kr-)f\:^ 
* 0 

< X,r ) (p 3 C x, r ) + w, (X, r ) r ) 

rdr 


“ N(A,) J[ W ' 




1 r dr 

J *? *7 


r Stx,) /E w ' 

( A <’^) <p ' 

(*4; 

) ♦ W* (*^) <p4\,^) 

1 r dr 

J*? *1 

&5/K 

■’i 

= n <t } - ) / L W ' 

Of, r ) <P, 


)+ w a (tj , r) <p a (\,^) 

J rdr 

C 5^X. 

e Ah/i; r 
= N(r,;/L W| 

7 o 

rV"> <P, 

( Ak, t 

)t W 2 (r f ,r) <P 2 ( x k4). 

j rdr 

&b / K 

1 

N(*p 

1 Eu 
1 E u 

* + v u P f 
i +v u J |_ 

o 

~ W 3 ( 

r) 



+ W*| 


) <P.(T K 

’ r L 

r dr 

rrj rrj 


c 

^'K 

_ e' 7 *'* 1 E.U 

W ( Ap EL l 

-v./f 
«+u J L 

_ w 3 1 

Irjj $3 ( ^<.7 r ) 



(4.29a) 
(4.29b) 
(4.29c) 
(4. 29d) 
(4.29e) 
(4.29f ) 
(4.29g) 
(4.29b) 

(4.291) 


+ c p*(r K , r)]^ ^ 


(4.29 j ) 



Integration of the functions G 2 j and G 3 ^ leads to 


\ 

2>i 


[2 v l t, (Yj) T t orj'rp -2Hi rtrp 


- %- 07 J09) ijr ^) + j;ctjo 


(4. 3°a) 


Vl P El J_ 

*V N C\/> EL l LU Euf 



(4.30b) 
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Discussion of Numerical Results . 

To demonstrate some of the aspects of the solution, 
selected numerical results are presented graphically in 
Figures 2 through 9* In the graphs which resulted from 
this numerical study, it is noted that the stresses are 
given in dimensionless form. Since A L is the average axial 
stress in the lower cylinder, the stresses are plotted as 
ratios CT Z /A L , T/A L , etc. The deflection, 6 , has been 
arbitrarily chosen as 


c l l\ h E l 1 ‘ 

Ei[ l Eu 7 ? J 


(4.30) 


The numerical study showed that the series, trun- 
cated after 35 terms, exhibited large oscillations for the 
stresses at the interface due to the discontinuous nature 
of the boundary conditions at the interface. The ampli- 
tudes were greatly reduced by using Feje'r's averaging 
method, as described by Lanczos [ 5 ]- on the first fifteen 
terms of the series. 

Figures 2 through 6 are for a configuration where 
^ J and the material constants and lengths are equal. 

The axial stresses as a function of r for problem 1 
is shown in Figure 2. The oscillations in the computed 
results at the interface are not unexpected. The sharp 
rise as the value r = 7? is approached is seen to have the 
same characteristic as equation (4.1) for the rigid flat- 
ended cylinder on an elastic half space. At z = 0.8 for 



the lower cylinders it is noted that a small outer portion 
is in tension. 

Figure 3 shows a schematic of the deformed grid for 
problem 2. The underlying characteristic is that at the 
interface, very small radial displacements occur compared 
with the rest of the cylinders. This is due to the re- 
straining action of the free outer annulus of the lower 
cylinder upper face. The portion of the grid shown dashed 
is not known exactly because of the lack of enough grid 
points. . The axial stresses at and near the interface is 
shown in Figure 4. The oscillations are more pronounced 
than in problem 1; otherwise, the characteristics are the 
same. Figure 5 shows the shear stresses near the interface 
and Figure 6 shows the radial and tangential stresses near 
the interface. 

A second configuration has been studied for pro- 
blem 2. It is an equal diameter case with E y = 0.2E L . 

For this case. Figure 7 shows the deformed grid. It is 
noted that, in view of the different E values, the vertical 
grid lines are not smooth through the interface. The 
slight curvature of the deformed interface is caused by 
the shear stresses transmitted through the interface. 

A limited study was then made to determine how the 
overall axial stiffness is affected by changing the ratio 
of cylinder radii. The solution from problem 2 was used 
for this study and results were plotted as dimensionless 
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ratios in order to be applicable to any radius, R L , for the 
lower cylinder. The value, "*) , is then the ratio, Ru/R l . 

Figure 8 represents a stepped cylinder configuration 
where the length of the lower cylinder was taken to be 
equal to its diameter. It shows A L as a function of "> 7 . 

The ordinate, A L , has been defined in this figure as 

A L = Fi/ir«ER L /(l + R l 2 Au) (4.31) 

It is noted that the curve begins and ends with a value of 
1.0 and exhibits a minimum value near R u /R l = O.32. When 
Ru/R L = 1.0, the cylinders are in a state of uniform axial 
stress. For this reason, the curve represents the error 
factor if the axial stiffness, F/'g , for a stepped cylinder 
is computed on the uniform axial stress assumption, as 
commonly done in strength of materials methods. It is seen 
that, for this configuration, the error can be as much as 
7 % (near Ru/R l = O.32). 

Figure 9 represents a stiff cylinder bonded to a 
cylinder having the characteristics of rubber. It shows 
A L /<5 as a function of Tj , The ordinate, A L /g, is seen to be 

A l /<5 = F 1/ 7r 5 E L R ^ (4.32) 

Two curves for two typical values of the ratio, &/R l , is 
given. For ^/R L =2.0, it is seen that the curve starts at 
zero, increases to a maximum near the value Ru/Rl = 0.95, 
and then drops slightly as the value Ru/R l = 1.0 is ap- 
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proached, That the curve exhibits a maximum is a paradox 
which still remains unexplained. For i/R L = 0.5» the 
curve has similar characteristics except that it is not 
certain if a maximum value is present. 
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t/a l 



Figure 5* Shear Stresses near the Interface, 

Problem 2, 7? = 0.5, 

h = £ = 1.0, V L = V u = 0.3, E L = E u s= 1.0 
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cr/A L 


Radial stress, tf r /A t 



Figure 6. Radial and Tangential Stresses near the Interface, 
Problem 2, 1) = 0.5, 

h m jt = 1.0, V L =Vu = 0.3, E L a* E y = 1.0 



Figure 7. Deformed 

h - i - : 



Figure 8. Axial Stiffness, Problem 2, Stepped Cylinder 
Configuration, V L = = 0.3» E L = E y , £ - h 



Figure 9« Axial Stiffness , Problem 2, Stiff Upper Cylinder 
Bonded to a Rubber Cylinder, V L = 0„5» 

Vy = 0.3> Eu = 100E t 




CHAPTER V 


SUMMARY AND CONCLUSIONS 

The solution to a two cylinder contact problem was 
achieved by utilizing the known solution of the finite 
circular cylinder and applying the appropriate boundary and 
continuity conditions. The interrelated Fourier coeffi- 
cients were expressed implicitly by an infinite system of 
linear algebraic equations. By truncation, an explicit 
solution of the Fourier coefficients was obtained. Two 
example problems were then solved. It was shown that a 
part of the uniform axial stress solution is not known until 
the Fourier coefficients were obtained. 

A typical cylinder configuration was chosen for a 
numerical study of each problem. It was noted that appre- 
ciable oscillations in the calculated stresses occurred 
at the contact interface of the cylinders. These dimin- 
ished rapidly a short distance away from the interface. 

The comparison of calculated results with those known for 
the rigid punch seem reasonable. For this reason, the 
exact nature of the stress singularity due to the discon- 
tinuous nature of the boundary conditions could not be de- 
termined. 

A definite determination of the validity of this 
solution is also precluded for this very reason. Since 



the equilibrium and compatability equations, Hooke's law, 
and the stress free boundary conditions on the curved sur- 
faces are satisfied implicitly from the eigenfunctions, 
the solution is valid if the rest of the boundary condi- 
tions are met. Those at the lower and upper rigid faces 
are met identically. At the interface, the boundary con- 
ditions are met in an averaged or integrated sense. The 
numerical results appear to provide useful information. 

To demonstrate the practical aspects of the solution 
a limited study was made to determine the effect of cylin- 
der radii ratio on the constants, A L and A L / 5 . It was 
noted how these constants can be used to obtain the true 
deflection rate of two cylinders bonded together. Some 
examples of the applications that this solution could be 
applied to ares shoulder bolts, integral pistons and rod, 
vibration and shock isolators, and concrete piers. 

The numerical results that were presented were ob- 
tained by using the double precision mode on a SDS Sigma 7 
computer. Seventeen pairs of eigenvalues were used. In- 
tegration of the matrix elements was obtained in closed 


form. 
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TABLE OF FUNCTIONS 


The <p functions used in this thesis are? 


<p, { r j>b) = W ~ V* W r j: (¥) 

•p.^fc) = [io-v)i;-j;(if)+ i(^) s r x(^f) 

%(>?. £) * Vitr,) r x(*f)- [ifl-vjr,- Xfv.o - r^T.c^ r^) 

<Psr ( Y ' , s) = [V * w + 1 s r°cv\ dfe) + Vjo-vi £ x(%f) 

- [zo-v) r,xfr;) + r, 1 jifr,)] 5. x,(3f) 

<P t (•»■/,£) = [zo-v) T,- X(r,) +• T ( z XCt,)J S. 

- f I - 2 -UJ) r/ J7Cn) 



k5 


The W functions ares 
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